Electronic transport through domain walls in ferromagnetic nanowires: 
Co-existence of adiabatic and non-adiabatic spin dynamics 
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We study the effect of a domain wall on the electronic transport in ferromagnetic quantum wires. 
Due to the transverse confinement, conduction channels arise. In the presence of a domain wall, spin 
up and spin down electrons in these channels become coupled. For very short domain walls or at high 
longitudinal kinetic energy, this coupling is weak, leads to very few spin flips, and a perturbative 
treatment is possible. For very long domain wall structures, the spin follows adiabatically the local 
magnetization orientation, suppressing the effect of the domain wall on the total transmission, but 
reversing the spin of the electrons. In the intermediate regime, we numerically investigate the spin- 
dependent transport behavior for different shapes of the domain wall. We find that the knowledge 
of the precise shape of the domain wall is not crucial for determining the qualitative behavior. 
For parameters appropriate for experiments, electrons with low longitudinal energy are transmitted 
adiabatically while the electrons at high longitudinal energy are essentially unaffected by the domain 
wall. Taking this co-existence of different regimes into account is important for the understanding 
of recent experiments. 

PACS numbers: 72.10.-d, 75.47.Jn, 73.23.-b, 75.75.-|-a 
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I. INTRODUCTION 

A new kind of electronic devices taking advantage of 
the electron spin have been developed during the last 
years. The influence of the spin on electronic transport 
attracts considerable interest since early experiments in 
multi-layered magnetic structures have shown that the 
resistance is considerably increased in the case of an anti- 
parallel magnetization of the layers, as compared to a 
parallel configuration^'^. This is at the base of the so- 
called Giant Magneto- Resistance (GMR), which is al- 
ready used in the read-heads of commercial high per- 
formance hard-disks. 

In magnetic configurations that arc obtained when 
one substitutes the non-magnetic spacer layer between 
the ferromagnets by domain walls'', the effect of a mag- 
netic domain wall on the electronic transport proper- 
ties has become a subject of great interest. In partic- 
ular, the effect of a single domain wall on the resis- 
tance of a ferromagnetic nanowire has been measured 
for electro-deposited cylindrical Co wires down to 35 nm 
in diametei— , and thin polycrystalline Co films having a 
thickness of 42 nm and a width down to 150 nm^. The re- 
sults of both experiments indicate that, besides negative 
contributions from the anisotropic magneto-resistance, 
the domain wall scattering yields a positive contribution 
to the resistance. 

The prospect of interesting technological applications 
of magneto-electronic devices exploiting the spin degree 
of freedom of the electrons together with its importance 
from the fundamental point of view has strongly mo- 
tivated theoretical studies of spin-dependent electronic 
transport. 



Many efforts have been made in order to explain 
the electronic transport, in particular the enhancement 
of the magneto-resistance in these spin-dependent elec- 
tronic devices. For example, a Boltzmann equation 
has been applied to study the resistance of multi- 
layered magnetic/non-magnetic structures when the 
spin-diffusion length is larger than the mean free patbS. 
While the electronic spin is expected to follow adiabati- 
cally a very slowly varying magnetization^, the deviations 
from this adiabatic behavior, which are due to the finite 
length of the domain wall, lead to a so-called mistracking 
of the spin, and result in a GMR-like enhancement of the 
magneto-resistance in multi-domain wall configurations^. 
Including spin-dependent scattering, the mistracking and 
the resulting magneto-resistance have been calculated for 
such a system, within the so-called two-band model«, 
which consists in a simplification of the complicated band 
structure of a ferromagnetic metal. 

An outstanding problem in magnetism is a fully con- 
sistent description of transport and thermal properties in 
terms of electronic states calculated from first principles. 
Despite remarkable progress in the past twenty years, 
such a description does not yet exist in a form suitable 
for predicting features such as domain wall structures 
in non-equilibrium situations. As such, it is reasonable 
to search for suitably simplified model descriptions that 
capture the essence of the important physics involved. 
In the present two-band model is useful for the 

study of how the geometrical characteristics of a mag- 
netic domain wall affect electron transport. We make a 
distinction between spatially extended electronic states 
that contribute strongly to conduction and more local- 
ized states that contribute strongly to the formation of 
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local magnetic moments. For the transition metals, this 
model assumes that exchange correlations between elec- 
trons in primarily d-like orbitals are largely responsible 
for the formation of magnetic moments leading to the mi- 
croscopic magnetization. The s-like orbitals contribute 
much less strongly to the magnetization and instead in- 
teract relatively weakly with the local moments via a 
contact interaction term. These types of s-d interaction 
models have proven very useful in the past for discussions 
of indirect exchange interactions in magnetic transition 
metal multilayers. 

In our work we therefore assume that the domain wall 
represents a stable magnetic state with an energy above 
the ferromagnet ground state. The exact shape and di- 
mensions of the wall are determined by exchange corre- 
lation energies and spin orbit interactions primarily af- 
fecting electrons associated with the magnetization of the 
wall. These interactions are small perturbations on the 
conduction electron states, and the interaction between 
conduction electrons and the domain wall is represented 
by a simple contact potential. In a single electron picture, 
the wall appears as a spatially varying spin dependent 
potential for the conduction electrons. The magnitude 
of the splitting between the spin up and spin down po- 
tentials is taken as a free parameter that is related to 
the exchange correlation energy of the electrons involved 
in forming the wall, and the contact potential describing 
interaction of conduction electrons with the effective po- 
tential associated with the wall structure. In the follow- 
ing we refer to this contact potential between conduction 
electrons and the magnetization as an 'exchange inter- 
action' although it is quite distinct from the exchange 
interaction used to parameterize the interactions leading 
to magnetic ordering. 

In thin ballistic quantum wires and narrow constric- 
tions or point contacts, the lateral confinement of the 
electronic wave- functions leads to the emergence of quan- 
tized transport channels. As a consequence, the conduc- 
tance is quantized and exhibits steps of /h as a function 
of the Fermi energji^iiSi. Nakanishi and Nakamuraii^ con- 
sidered the conductance of very narrow quantum wires 
including the effect of a domain wall. A perturbative 
approach allowed them to study the effect of a very 
short domain wall on the conductance steps. Imamura 
and coUaboratorsi^ calculated the conductance of a point 
contact connecting two regions of a ferromagnet having 
parallel or anti-parallel magnetization directions. Within 
an s — c? two-band model, they numerically obtained a 
non-monotonic dependence of the domain wall contribu- 
tion to the resistance on the width of the point contact. 

There have been attempts to compare different ap- 
proaches to the calculation of the domain wall magneto- 
resistance (DWMR), which point towards the importance 
of including more realistic band structures. In the bal- 
listic case, van Hoof and collaborator»i2i calculated the 
DWMR for an adiabatic model where the magnetization 
direction changes very slowly along the wire, using an 
extension of the standard band structure calculation to 



include an infinite spin spiral, as well as a "linear" model, 
where the magnetization turns in a finite region at a con- 
stant rate. These two approaches yield a much larger 
effect than a two-band model, where corrections with re- 
spect to an infinitely long domain wall are calculated. On 
the other hand, first principle calculations for the case of 
abrupt magnetization interfaces yield a DWMR which is 
orders of magnitude larger than for the other models that 
take into account realistic domain wall lengths. 

It then seems necessary to develop more accurate treat- 
ments within the two-band model in order to under- 
stand the crossover from the abrupt domain wall sit- 
uation to the adiabatic regime. Since we are working 
with nanowires, obvious transverse quantization effects 
appear, which are more easily tractable within a two- 
band model. Moreover, the two-band model allows to 
easily obtain the transmission coefficients with and with- 
out spin flip, making it possible to study the mistrack- 
ing effect in finite length domain walls and its GMR-like 
consequences. Finally, given the typical experimental pa- 
rameters, it would be important to go beyond the ballis- 
tic limit. Taking into account disorder within a two-band 
model seems much more doable than in the framework 
of a band structure calculation. 

In this paper, we study the effect of a domain wall 
on the conductance of a nanowire within the two-band 
model, comparing different shapes and sizes of the do- 
main wall. Focusing on the contribution of the spin- 
dependent scattering of the domain wall, we do not con- 
sider material-dependent contributions to the resistance 
like the anisotropic magneto-resistance. 

After presenting our model in sectional we consider 
the perturbative regime of weak spin coupling induced 
by the domain wall in section UTTI From a comparison of 
different domain wall shapes, we shall extract the rele- 
vant parameters governing the transmission through the 
domain wall with and without spin-fiip processes. While 
this applies to short domain walls, in section lTvl we study 
the transport in the general case of a strong spin coupling 
induced by the domain wall, which allows to treat domain 
walls of arbitrary length. We present the deviations from 
the full transmission with spin rotation following the lo- 
cal magnetic structure expected for infinitely long do- 
main walls, which is due to finite domain wall length, 
and treat the most interesting crossover regime. This 
case is relevant since typical experiments^ are far from 
the thin wall regime, but not really in the adiabatic limit 
in which the length of the domain walls is large. 



II. MODEL 

We consider a wire along the z axis with a domain wall, 
and choose the origin {z = 0) in the middle of the wall 
(see Fig. As explained above, we work within the 
two-band model, where the d electrons are responsible 
for the magnetization and the current is carried by the s 
electrons. Therefore, we write for the latter an effective 
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FIG. 1: Sketch of a ferromagnetic quantum wire containing 
a domain wall (grey region), for the example geometry of a 
square cross-section. The arrows indicate the magnetization 
directions far from the domain wall, for the case of a Neel 
wall. 



Hamiltonian 

where A is the spin splitting of the s electrons due to the 
exchange coupUng with the d electrons and a is the vec- 
tor of the Pauli matrices. The unit vector / represents 
the direction of the local magnetization. Its functional 
dependence describes the shape of the domain wall. The 
lateral confinement present in a nanowire may have a con- 
siderable influence on this shape, leading to domain walls 
which are altered as compared to the case of bulk do- 
main walla^ii^ii^iiSiii. In addition, a spin-polarized cur- 
rent through the domain wall creates a torque which can 
alter its shapeiS.. Working in the linear response regime 
of low current, we do not need to take into account this 
back-action of the conduction electrons on the magnetic 
structure. 

Assuming that the magnetization only depends on z, 
we will therefore be interested in comparing different 
functional forms of the kind / = {fx{z),0, fz{z)}. This 
choice does not imply a loss of generality and corresponds 
to a so-called "Neel wall", where the magnetization is 
parallel to the wire axis in the leads far from the domain 
wall (arrows in Fig. and turns inside the x-z plane 
parallel to this wire axis when going through the wall. 

The assumption f{r) = f{z) allows us to separate the 
transverse and longitudinal parts of the Hamiltonian 
The transverse quantization gives rise to transport chan- 
nels with quantum numbers ni and 712, and an energy 
En-i,n2- The density of channels p = 2T:mA/h^ (where A 
is the cross-section of the transverse area of the wire) is 
equal to that of a two-dimensional system, and thus inde- 
pendent of the energy and the shape of the cross-section. 
For the example of a wire with a square cross-section of 
side w, we have 




Far away from the domain wall, the orbital parts of 
the eigenstates are products of transverse channels and 



longitudinal plane waves. Since lim^^±oo /^(z) = ±1, 
the associated eigenenergies for spin up are 

i^T=^n,,„. + ^±y, (3) 

while for spin down we have 

The domain wall leads to the scattering of these states 
and the conductance (in units of e^/h) is given by the 
Landauer formula 

5= E E^"i"".(^f)' (5) 

ni ,712 (7,{7' 

where the sum is done over the occupied channels. T^fn2 
is the transmission coefficient in the channel (rti, 712), for 
scattering of electrons with spin a into spin cr'. Such a 
coefficient only depends on the longitudinal energy 

e = Ep — Eni.n2 (6) 

as 

nfUEF)^T'^-^'{e). (7) 

Therefore, for each channel (ni,n2), we are left with an 
effective one-dimensional problem at energy e. 

In order to determine the transmission probability 
T'^''^ (e), we write the spinor wave- function in the up- 
down basis (with fixed, z-indcpendcnt spin orientations) 
as 

IV'(z)) =0TWk,T) + 0iWk,i), (8) 

where |z, ]) has to be interpreted as the tensor product 
of the position eigenvector \z) and the spin up state | t 
). The Schrodinger equation with the Hamiltonian 
leads to a system of coupled differential equations for the 
components 4>^{z) and 4>i{z): 

cP ^ 2m f A \ , 2m A , 

While the term containing plays the role of a spin- 
dependent potential, the transverse component fx of the 
wall profile is responsible for the coupling between the 
spinor components (j)^ and . The scattering solutions 
of Eq. are then needed to calculate the transmission 
coefficients, and therewith the conductance through the 
domain wall. 

For the extreme case of an abrupt domain wall, when 
fz has a jump from —1 to 1, the right-hand-side of Eqs. 
@ vanishes, and spin up and spin down electrons remain 
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uncoupled. The only effect of the discontinuity of fz is 
a spin-dependent potential step of height ±A for spin 
up/down electrons. Incoming spin up electrons having 
longitudinal energy e < A/2 cannot overcome this step 
and are reflected with probability one. Since the den- 
sity of conduction channels is independent of the energy 
for wires having a two-dimensional cross-section, this 
mechanism blocks a fraction A/2Ep of the conduction 
channelsi^, all of which exhibit perfect transmission in 
the absence of the domain wall. If one neglects the effect 
of the potential step on electrons having higher longitu- 
dinal energy (e > A/2), this channel blocking mechanism 
leads to a relative change in conductance 



due to the presence of the domain wall. Taking into ac- 
count the spin conserving reflections for e > A/2 leads^^ 
(in the limit Ep ^ A) to an increase of the effect by a 
factor 4/3. 

The precise shape of the actual domain wall present 
in an experimental measurement (which is very difficult 
to know) would in principle be needed to determine the 
scattering states. In addition, it is not possible to find an 
analytical solution of Eqs. for arbitrary domain walls. 
This is why we introduce various models of a domain 
wall, and approximate analytical, as well as numerical 
calculations. 

In the bulk, when the magnetization always remains 
parallel to the domain wall, we have the so-called Bloch 
walls, whose shape was originally calculated by minimiz- 
ing the total free energy in the thermodynamic limilj^i. 
In this case, /(r) is given by 

/(z) = {tanh(|),sech(|),0} , (11) 

where A is the length scale of the domain wall. The lat- 
eral confinement present in a nanowire will certainly al- 
ter the previous functional form of /. Moreover, an easy 
magnetization axis in the direction of the nanowire will 
result in a Neel wall, modified by the transverse confine- 
ment. Such effects have been recently discussed in the 
literatureiSiiSiii. 

The variety of possible domain wall structures moti- 
vates us to consider different domain wall profiles: linear, 
trigonometric and extended (defined below), in order to 
determine the influence of the domain wall shape on the 
conductance of the wire. As we will see below, while most 
of the effects are not very sensitive to the details of the 
wall, the signature of the particular domain wall appears 
in some regimes. 

A possible starting point is to assume that the Neel- 
like confined domain wall has components with the same 
functional form as in Eq. (|ll|l . In this case we will con- 
sider the "extended" domain wall defined by the magne- 
tization direction 

/l-)(^) = {sech(|),0,tanh(^)} . (12) 



As compared with the situation of a Bloch wall, this leads 
to a permutation of the spatial variables in Eqs. ©, and 
does not change the results for the transmission coeffi- 
cients. This is why in the limit of high electron energy 
we can compare our results with the ones of Cabrera and 
Falico\«2^, who considered Bloch domain walls. 

In the case of weak coupling between the spin up and 
down states described by Eq. Q and short domain walls 
it is reasonable to approximate fz in the wall profile with 
a linear function of position 

P"'Hz) = I {VWW, 0, z/x] , for \z\ < A 
\ {0,0, sgn(z) } , for \z\ > A. 

(13) 

The semi-circle form of the coupling term in this "lin- 
ear" wall ensures that |/(z)P = 1. This is an important 
difference as compared to our previous work^^, where the 
conservation of the absolute value of the magnetization 
was not respected. The above constraint has only quan- 
titative consequences in the short wall limit which play a 
role when comparing different wall profiles, but becomes 
crucial for longer domain walls in the adiabatic regime. 

For an arbitrary domain wall, the extension of the stan- 
dard recursive Green function methodSiSi to take into 
account the spin degree of freedom (in a tight binding 
setup) allows us to calculate the transmission and reflec- 
tion coefficients Tn, T||, and R^i- 

The case of a "trigonometric" domain wall 

/^*")rz^ - / 0' sinff } , for \z\ < A . . 

^ " \ { , , sgn(z) } , for \z\ > A ^ > 

admits an exact solution for the wave-function inside the 
domain wall^^i which has recently been used to calcu- 
late the torque that is due to a spin-polarized currenti^. 
We use the exact solution to determine the scattering 
properties of the domain wall. Details are presented in 
Appendix 1X1 

By comparing with this exact analytic solution we 
checked the accuracy of the numerical method, as well 
as the absence of lattice effects for the parameters that 
we work with. 



III. WEAK COUPLING 

The problem can be treated at the analytical level 
when the differential equations ^ are only weakly cou- 
pled. This is the case when the domain wall (in which 
the spin-flip terms fx are non-zero) is very short, or when 
the longitudinal energy of the electron is very high such 
that the transmitted electrons spend only a short time 
inside the domain wall region. 

Then, we treat the system of coupled differential equa- 
tions lO iteratively, considering the spin-flip terms on the 
right-hand-side as a perturbation. Within this approach, 
the solution to the homogeneous differential equation 
(|9a|l (in which the spin-flip terms induced by f^ are set to 
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zero) is injected in the spin-flip term of the second differ- 
ential equation H9b|l. This method was used in Ref. 
where a mechanism of channel blocking by a domain wall 
was proposed as a source of resistance in short ferromag- 
netic quantum wires. The starting point of this approach, 
which we present here for the example of a linear domain 
wall as described by H13|) . is an incoming majority (spin- 
up) electron from the left (p^, and (j>^ = 0. Outside the 
domain wall region, (pf reads 



■ikz 



with the wave-numbers 
k = 



ik' z 



for 
for 



z < -A 
z> X, 



1 2m 



A 



(15a) 
(15b) 



(16a) 
(16b) 



For — A < z < A, the homogeneous solution of Eq. (|9a|l is 

2e 



(t)f{z) = aAi 
-t- /3Bi 



2/3 _ 



2/3 



Z 

A ■ A 

2e z 
A + A 



(17) 



with the usual Airy functions Ai and Bi, and the dimen- 
sionless parameter p defined by 



1/2 / aY^' 

A = ( TTT^ I fcpA . 



\2EfJ 



(18) 



The coefficients a and /3, as well as if| and r||, are ob- 
tained from the matching of the expressions (|15|l and 
()17|l at z = ±A. For e < A/2, we have imaginary k'. The 
transmission without spin- flip is zero and \r-\i\ — 1. For 
e > A/2, the transmission without spin-flip is finite. 

The first-order correction from the spin-conserving 
scattering is then obtained by injecting (j)^ as the inho- 
mogeneous term in the differential equation (|9b|l for (j)^. 
Since we do not have incoming spin-down electrons, out- 
side the domain wall region we take the outgoing plane 
waves 



— ik'z 



ikz 



for 
for 



z < -A 
z > A . 



(19a) 
(19b) 



For — A < z < A the general solution can be written as 
a linear combination of Airy functions plus a particular 
solution (/)P(z) 



ai Ai 
/3i Bi 



+ 0P(z). (20) 
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FIG. 2: (a) Perturbative results for T-^ from Eqs. II22II and 
12311 (solid lines) for three values of p, compared with the 
corresponding full numerical results (circles, squares and dia- 
monds) for the linear domain wall, (b) The total transmission 
T and Tn (diamonds and triangles, respectively) for p — 0.20. 



Spin-flip processes are now included in the description, 
and electrons undergoing a spin-flip can be transmitted 
even for energies e < A/2. The parameters ai and /3i 
and the coefRcients t^ ^ and r^i are determined from the 
matching of and IpH) . 



A. Long wavelength limit 

In the limit of a short wall, when the wavelength of the 
incoming electron is much longer than the domain wall, 
the linear approximation of the Airy functions allows us 
to write for e > A/2 



ItTTi' 



4kk' 



(k + k'Y + A{\kk'Y 
{k - k'f + ^{Xkk'f 

{k + k'Y + AiXkk'Y 



(21a) 
(21b) 
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Obviously, for A — > we recover the well-known results 
for a step potential^^. 

For e < A/2, the transmission probability with spin- 
flip is given by 



n 



\tn\' 



cV(i + | 



(22) 



For e > A/2, 



Jti = 40 p 5— 

((V^TT + V^^)2 + V(x2 - 1))' 

(23) 

is obtained, with x = 2e/A and the prefactor C defined 
by 



C 



(24) 



In Fig. m (a) we show from Eqs. and as 
a function of 2e/A, together with numerical calculations 
for three different values of p. We can see an excellent 
agreement for the smallest values of p. When the value 
of p is increased, deviations appear first at energies close 
to A/2. These features are consistent with the fact that 
the linear approximation of the Airy functions is justified 
in the small p limit, and becomes increasingly better for 
small energies. 

The total transmission T ^ Ti^i + (Fig. H (b)) 
is dominated by the large transmission without spin-flip 
for £ > A/2. This feature justifies the "channel block- 
ing" picture proposed for short domain walls in Ref . |l9j | , 
where the presence of the wall suppresses almost com- 
pletely the transmission at energies e < A/2 {T^i is 
only a small correction for e < A/2 and negligible for 
e ^ A/2). Neglecting what happens for e > A/2, we 
have Sg/g = (-1 + {Cp)^)A/2Ef. 

In the limit of short domain walls, we have found that 
the dependence of the transmission coefficients on the 
shape of the wall is only through the integral over 
which enters in the prefactor C. For the linear, trigono- 
metric and extended domain walls, C takes the values 
7r/2, A/tt and tt, respectively. Such a scaling is shown in 
Fig. 131 where the transmissions T^^ divided by the cou- 
pling strength {CpY for the different domain wall shapes 
coincide for all energies, except those close to A/2. 



B. Short wavelength limit 

The perturbative approach is not only applicable for 
short walls and low energies (as in section UlI A|l . but for 
general domain wall parameters as well, provided that 
e ^ {AX/h)^ /m. That is, when the time that the elec- 
tron spends inside the domain wall is much shorter than 
the spin precession period, and therefore the spin-flips are 
very unlikely. In this limit, the WKB approximation of 
the scattering wave-functions for the linear domain wall 
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FIG. 3: The transmission divided by the coupling 

strength (Cp)^ for the linear, trigonometric and extended do- 
main walls (solid line, up and down triangles, respectively). 
The inset shows the results for T|| for the same value of 
p = 0.09, before dividing by the corresponding value of (Cp)^. 
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FIG. 4: The coefficients Tf^ and J?|| as a function of Cp for 
the three different domain wall shapes, at e = 0.95A/2. An 
oscillatory behavior for the linear and trigonometric walls is 
found as a consequence of the edges at the connection to the 
leads. 



model (Eq. yields the reflection and transmission 
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FIG. 5: Transmission and reflection coeflicients in the inter- 
mediate regime (p = 1) for the extended domain wall geom- 
etry 11211 . A large transmission with spin-flip (diamonds) is 
found for energies e < A/2 in this regime, where the transport 
is adiabatic. 



coefficients 



Ay sin2(2fcA) 



Ti 



TT 



SI 

1 - - 



sin2(2fcA) 



Ti 



(25) 

(26) 

(27) 
(28) 



Thus, for energies e ^ A all scattering coefficients, ex- 
cept the transmission without spin- flip, are very small. 
Therefore, in first approximation we can neglect the effect 
of the domain wall for electrons with high longitudinal en- 
ergies. The conductance associated with the domain wall 
is then determined by the low-energy electronsiS.. The 
algebraic decay in A/e is less pronounced than the expo- 
nential suppression obtained by Cabrera and Falico\iS^. 
Such a difference arises from the sharp edges at z = ±A in 
the linear domain wall model we used for this calculation. 



IV. STRONG COUPLING 

If we are interested in energies e ~ A/2 and not neces- 
sarily short walls, the previous picture has to be modified. 
The linear approximation and the perturbative treatment 
(involving only one spin-flip) in the wall region are no 
longer justified. Beyond the perturbative regime, the de- 
tailed shape of the domain wall might become relevant. 

In Fig.^ we present and as a function of the 
coupling strength for different domain wall shapes and 



FIG. 6: The relative change in conductance 5g/g caused by 
the presence of a domain wall (extended shape), for = 2 A. 



an energy of the order of A/2. As discussed in the previ- 
ous section, it is for these energies that a dependence on 
the detailed shape of the domain wall appears first when 
departing from the weak coupling limit (Cp <C 1). We 
can see from Fig. 0]that the transmissions (reflections) 
for the different domain wall shapes coincide for small 
values of Cp. Even for stronger couplings, the different 
models do not show very important differences in their 
behaviors. The only apparent difference are oscillations 
of the transmission coefficients as a function of p, which 
occur at intermediate p for linear and trigonometric do- 
main walls. The origin of these oscillations is due to the 
edges of the domain wall region leading to Fabry-Perot 
like interferences. For a smooth domain wall structure 
such as the extended domain wall, the oscillations are 
absent. On the other hand, and as expected, T^^l 1 for 
all shapes in the limit of large p. 

It is the limit of infinite domain wall length where 
the spin follows adiabatically the orientation of the local 
magnetization^, corresponding to a rotation from spin 
up to spin down in the external basis of fixed spin orien- 
tations. Electrons are transmitted with probability one 
through the wall, therefore T-\^i — 1 and T|| — = 
R-\l = 0. In this limit the detailed shape of a domain 
wall, having slow spatial spin rotation, is irrelevant. The 
adjustment of the spin to the direction of the local mag- 
netization requires an infinite number of spin-flips (in the 
fixed basis), and obviously cannot be described by tak- 
ing into account a small number of spin-flips as in the 
perturbative approach used for short domain walls. The 
condition for the local adjustment is that the Larmor 
precession of the spin around the local magnetization is 
much faster than the rotation of the local magnetization 
viewed by the traveling electron^. This condition of adi- 
abaticity translates into A ^ {h/Xj^/TJm. 

We then see that the adiabatic condition strongly de- 
pends on the longitudinal kinetic energy of the electrons. 
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In a quantum wire, at a finite value of A, electrons with 
low longitudinal velocity are essentially adiabatic, while 
the channels with low transverse quantum numbers can 
be highly non-adiabatic. In calculating the conductance 
of a ferromagnetic quantum wire, we have to take into ac- 
count the co-existence of adiabatic (low longitudinal en- 
ergy) and non-adiabatic (high longitudinal energy) elec- 
trons. It then seems important to work out the crossover 
between the short wall and adiabatic limits, for different 
shapes of the domain wall. 

For an intermediate value of p in the case of an ex- 
tended domain wall, Fig. [S] shows that the behavior for 
e < A/2 is radically different from the weak coupling 
case of section urn The weak coupling result of Eq. H22|) 
is only valid at extremely low energies, and the Ti^i ap- 
proaches one (adiabatic behavior) for longitudinal ener- 
gies considerably lower than the step height. Above A/2, 
T|j decreases monotonously with energy, returning to the 
weak coupling regime in the limit of large e. At the same 
time, increases towards one and T remains very close 
to perfect transmission for all energies. 

Thus, almost all of the electrons with energy e > A/2 
are transmitted. However, while the spin of the trans- 
mitted electrons is changed by the domain wall for low 
e and high p, the spin of electrons having high e in do- 
main walls of low p remains unaffected by the wall (see 
also Fig. [21(a)). Therefore, in calculating the effect of 
the domain wall on the quantum conductance the modes 
with longitudinal energies in the interval (—A/2, A/2) 
are most relevant^®. 

The conductance for the ideal ballistic case given in Eq. 
is obtained by summing over all conductance chan- 
nels. Fig. 1^ shows an example of the resulting behavior 
for the difference in conductance between the cases with- 
out and with domain wall (normalized to the conductance 
without domain wall), as a function of the domain wall 
parameter p. We can see that the channel blocking ef- 
fect due to the presence of the domain wall is rapidly 
suppressed upon increase of the coupling. Similar results 
have recently been obtained using a different numerical 
approaches. 



main wall is not relevant, and the transmission coeffi- 
cients scale with the coupling strength. For a very long 
domain wall the spin of the electrons follows adiabati- 
cally the local effective magnetization and the conduc- 
tance is unaffected by the domain wall, independently of 
its shape. 

The intermediate coupling regime is most relevant for 
the domain walls that can be investigated experimen- 
tally. We have shown that the degree of adiabaticity of 
electrons at the Fermi energy strongly depends on their 
longitudinal kinetic energy. While the spin of electrons 
with low longitudinal energy essentially behaves adiabat- 
ically, the spin of electrons with high longitudinal en- 
ergy is practically unaffected by the domain wall. The 
crossover between these two behaviors, as a function of 
the longitudinal energy of the electrons, has to be taken 
into account in calculating the conductance of the quan- 
tum wire. 

Our analysis has been based on coherent scattering 
at the domain wall, which is connected to perfect leads. 
However, in realistic situations, the domain wall is not 
connected to scattering-free regions. The imperfections 
and impurities at both sides of the wall give rise to elas- 
tic scattering, which may be different for the two spin 
directions of the electrons. Though these coherent ef- 
fects can in principle be taken into account in a scatter- 
ing approach, such a coherent picture is not sufficient for 
wires which are longer than the phase coherence length. 
Since this is the case in typical experiments, we need 
in addition to take into account inelastic processes (like 
electron-phonon or spin-magnon scattering). The length 
of the leads over which the spin of the electrons is con- 
served can then be described phenomenologically by clas- 
sical spin-dependent resistorsSi. In this situation, the 
important part of the electrons which do not undergo 
spin-flip processes leads to an increase of the resistance 
due to the GMR mechanism. This picture is likely to be 
representative of the experimental situations . More 
experimental and theoretical work concerning the various 
relaxation rates will be necessary to establish a complete 
quantitative understanding of the phenomenon. 



V. SUMMARY AND CONCLUSIONS 

The effect of a single domain wall on the electronic 
transport in a ferromagnetic nanowire has been studied 
systematically in various parameter regimes. The domain 
wall leads to a coupling of spin up and spin down elec- 
trons in the conduction channels, which is proportional 
to the exchange energy of the conduction electrons and 
the length of the wall. 

For an abrupt domain wall the step in the effective po- 
tential felt by the conduction electrons blocks the trans- 
mission of channels with low longitudinal energy. In the 
weak coupling limit, a perturbative approach is possible, 
leading to the lowest order correction to perfect chan- 
nel blocking. In this case, the detailed shape of the do- 
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APPENDIX A: EXACT SOLUTION FOR THE 
TRIGONOMETRIC DOMAIN WALL 



such that the solutions (C|,Cx) and 2me/h^ are the 
eigenvectors and eigenvalues, respectively, of the matrix 



A particularly instructive case is that of a trigonomet- 
ric domain wall (Eq. dJ) since an exact solution for 
the wave function inside the domain wall region can be 
obtained^^. Here we extend this approach to a scatter- 
ing situation by matching the inner solutions with plane 
waves, which allows us to calculate the transmission and 
reflection amplitudes. 

In addition to the external spin basis {|z, t), \z, J,)}, it 
is useful to introduce a local spin basis {|z, f'"), \z, i^)}, 
which corresponds to spin orientations parallel and anti- 
parallel to the (rotating) local magnetization direction 
/(z), leading to 



where the spin rotation matrix is given by 



cos(az -I- 7r/4) sin(az + 7r/4) 
- sin(az + 7r/4) cos(az 4- 7r/4) 



(Al) 



(A2) 



with a = 7r/4A. For z = — A, i? is simply the identity 
matrix (the local rotating basis coincides with the fixed 
one), and putting z = A corresponds to exchanging the 
spin directions between the local and fixed bases. 
Inserting the spinor 



iv;(z)) = 0;^(z)|z,f)-f 4(^)1^^ 



(A3) 



into the Schrodinger equation corresponding to the 
Hamiltonian we obtain 



dz2 



2m 



^ 2a 



dz 



(A4a) 



which is, in fact, Eq. (jSJ) expressed in the local spin basis, 
for the case of a trigonometric domain wall. This system 
of coupled differential equations can be reduced to a 2 x 2 
eigenvalue problem with the ansatz 



b]{z) 
h\{z) 



exp(ifcz) 



(A5) 



M 



—2ika 



2ika 

P + 0? - ra/\/h? 



(A6) 



The secular equation of M leads to the dispersion rela- 
tions 




4a2fc2 4_ ^ ( ^ 



(A7) 



that is, the eigenenergies of the infinite spin spiraP'^. For 
a closed spiral, the periodic boundary conditions would 
leadi to quantized values of k. However, we are interested 
in a scattering problem, where the region in which the 
magnetization turns is connected to homogeneous ferro- 
magnetic leads. We therefore express the general solution 
((/)| , 0| ) for a given energy as a linear combination of the 
four corresponding eigenstates of the spiral, and use the 
matching conditions between the domain wall region and 
the perfect leads at z = ±A. Taking into account the ro- 
tation of the local basis, and using the expressions given 
in Eqs. (|15|l and H19|l for the wave- function outside the 
wall region, we get 



ife'A 



ikX 



0L(A) 



ik (e 



-ikX ^ AkX^ 



dz 



0L(_A) 



dz 



a4>\{-X) 
a0L(A) 



ik' X 



-ik'r^ie 



_d_ 

dl 
d 



These conditions allow us to extract the amplitudes , 
r|-f, and r||, as well as the precise form of the wave- 
function inside the domain wall. The resulting transmis- 
sion coefficients are presented and discussed in Sec. IIVI 



^ M.N. Baibich, J.M. Broto, A. Pert, F. Nguyen Van Dau, F. 

Petroff, P. Eitenne, G. Creuzet, A. Friederich, J. Chazelas, 

Phys. Rev. Lett. 61, 2472 (1988). 
^ G. Binasch, P. Griinberg, F. Saurenbach, W. Zinn, Phys. 

Rev. B 39, 4828 (1989). 
^ J.F. Gregg, W. Allen, K. Ounadjela, M. Viret, M. Hehn, 

S.M. Thompson, J.M.D. Coey, Phys. Rev. Lett. 77, 1580 



(1996); M. Viret, D. Vignoles, D. Cole, J.M.D. Coey, W. 
Allen, D.S. Daniel, J.F. Gregg, Phys. Rev. B 53, 8464 
(1996). 

* U. Ebels, A. Radulescu, Y. Henry, L. Piraux, K. Ounad- 
jela, Phys. Rev. Lett. 84, 983 (2000). 

^ G. Dumpich, T.P. Krome, B. Hausmanns, J. Magn. Magn. 
Mater. 248, 241 (2002). 



10 



^ T. Valet, A. Fert, Phys. Rev. B 48, 7099 (1993). 

A. Stem, Phys. Rev. Lett. 68, 1022 (1992). 
* P.M. Levy, S. Zhang, Phys. Rev. Lett. 79, 5110 (1997). 
® B.J. van Wees, H. van Houten, C.W.J. Beenakker, J.G. 

Wilhamson, L.P. Kouwenhoven, D. van der Marel, C.T. 

Foxon, Phys. Rev. Lett. 60, 848 (1988). 
^° D.A. Wharam, T.J. Thornton, R. Newbury, M. Pepper, 

H. Ahmed, J.E.F. Frost, D.G. Hasko, D.C. Peacock, D.A. 

Ritchie, G.A.C. Jones, J. Phys. C: Sohd State Phys. 21, 

L209 (1988). 

" K. Nakanishi, Y.O. Nakamura, Phys. Rev. B 61, 11278 
(2000). 

H. Imamura, N. Kobayashi, S. Takahashi, S. Maekawa, 
Phys. Rev. Lett. 84, 1003 (2000). 

J.B.A.N. van Hoof, K.M. Schep, A. Brataas, G.E.W. 
Bauer, P.J. KeUy, Phys. Rev. B 59, 138 (1999). 
P. Bruno, Phys. Rev. Lett. 83, 2425 (1999). 

I. L. Prejbeanu, L.D. Buda, U. Ebels, K. Ounadjela, AppL 
Phys. Lett. 77, 3066 (2000). 

1® I.L. Prejbeanu, M. Viret, L.D. Buda, U. Ebels, K. Ounad- 
jela, J. Magn. Magn. Mater. 240, 27 (2002). 

^"^ B. Hausmanns, T.P. Krome, G. Dumpich, E.F. Wasser- 
mann, D. Hinzke, U. Nowak, K.D. Usadel, J. Magn. Magn. 



Mater. 240, 297 (2002). 

X. Waintal, M. Viret, arXiv:cond-mat/03 01293l 
D. Weinmann, R.L. Stamps, R.A. Jalabert in Electronic 
correlations: from meso- to nano-physics, ed. by T. Mar- 
tin, G. Montambaux, J. Tran Thanh Van (EDP Sciences 
2001). 

P.E. Falloon, R.L. Stamps, unpublished. 
L.D. Landau, E.M. Lifshitz, Course of Theoretical Physics, 
Vol. 8, 1 Ed. Pergamon Press (1960). 
22 G.G. Cabrera, L.M. Falicov, Phys. Stat. Sol. (b), 61, 539 
(1974); ibid 62, 217 (1974). 

P.A. Lee, D.S. Fisher, Phys. Rev. Lett. 47, 882 (1981). 
2" H.M. Pastawski, E. Medina, Rev. Mex. Fis. 47 SI, 

1 (2001), and references therein, also available as 
arXiv:cond-mat/0103219 

25 A. Brataas, G. Tatara, G.E.W. Bauer, Phys. Rev. B 60, 
3406 (1999). 

2® E. Merzbacher, Quantum Mechanics, John Wiley & Sons, 

2 Ed. (1970). 

2^^ Y. Imry, Introduction to Mesoscopic Physics, Oxford Uni- 
versity Press (New York 1997). 



